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—  This  report  describes  a  computer  program  and  its  operational 

details  which  transforms  stochastic  data  and  evaluates  probability 
distributions.  This  methodology  was  used  in  recent  papers  by  Snyder 
and  Thcmas,  1986;  and  Thcmas  and  Snyder,  1986a,  1986b. 
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INTRCDUCTION 

A  data  transformation  method  was  developed  and  tested  which 
appears  to  be  useful  with  a  wide  variety  of  stochastic  data  over  a 
wide  range  of  sample  characteristics.  The  transform  is  essentially 
exponential  and  allows  reflection  of  one-side  infinite  data  into  a 
finite  domain  with  proper  boundary  conditions.  The  mathematical 
transform  function  is  scaled  and  shaped  by  the  sample  mean,  standard 
deviation,  and  coefficient  of  skew.  The  method  provides  a  means  for 
visual  display  of  data  with  slight  curvilinearity  for  simple  samples, 
or  for  more  carp  lex  samples  exhibiting  such  characteristics  as 
bi-modal  tendency,  positive  or  negative  skew,  and  presence  or  absence 
of  outliers. 

In  an  earlier  research  report,  Thomas  and  Snyder  (1984c)  reported 
a  mathematically  form-free  method  for  performing  probablistic  analysis 
and  synthesis  of  experimental  data  without  assuming  any  particular 
conventional  distribution  function.  They  had  successfully  used  the 
methodology  in  studies  reported  in  three  papers  (Snyder  and  Thcmas, 
1983;  Thcmas  and  Snyder,  1984a,  1984b) .  Additional  background 
material  on  sliding  polynomials  which  is  used  in  the  related  smoothing 
operation  may  be  found  in  papers  by  Snyder  (1976  and  1980) . 

We  stated  in  the  earlier  report  (Thcmas  and  Snyder,  1984c)  that 
the  transform  would  be  gradually  standardized  through  additional 
research  and  experience.  The  need  for  improvement  of  the  transform 
became  very  clear  as  we  tried  to  develop  seasonally  continuous 
probability  distributions  (Snyder  and  Thcmas,  1986;  and  Thcmas  and 
Snyder,  1986b) .  The  old  version  of  the  transform  proved  inadequate  as 
we  tried  to  bring  into  alignment  the  samples  from  different  months  of 
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the  year  to  provide  nearly  homogeneous  data  that  could  be  smoothed  by 
a  mathematical  surface.  The  experience  and  success  of  the  improved 
transform  are  reported  in  the  following  three  papers:  Snyder  and 
Thomas,  1986;  and  Thomas  and  Snyder,  1986a,  1986b. 

The  description  and  computer  program  of  the  improved  transform 
are  provided  in  this  report  for  those  that  are  interested. 


4 


TRANSFORM  STRUCTURE 


The  mathematical  transform  consists  of  two  steps.  First,  the 
historical  variate,  h,  is  reduced  to  a  standardized  variate,  h'. 
Secondly,  h'  is  converted  to  an  abstract  variate,  v,  using  an 
ogee-shaped  curve  composed  of  two  limbs  of  the  carrmon  descending 
exponential,  joined  at  a  common  point,  C,  as  shewn  in  Fig.  1.  These 
two  limbs  are  both  shifted  and  shaped  with  the  skewness  of  the  sample. 
The  addition  of  the  skewness  parameter  improved  the  transform  and 
contributed  significantly  to  its  enhanced  capabilities.  Changes  in 
the  transform  for  positive  skew  are  shown  in  Fig.  1  while  changes  with 
negative  skew  are  idential  but  reversed. 

Positive  Skew.  The  v^  and  v^  limbs  of  the  transform  are  given  in 
equation  [1] . 

v.  =  VB  -  (VB  -  C)  exp  (f  h ' ) ,  v  >  C 
1  [1] 
v2  =  C  exp  (~d  h') ,  v  <  C 

where  C  is  the  common  point  of  limbs  in  v-scale,  VB  is  the  asymptotic 

right  boundary  of  the  v^  limb,  f  and  d  are  shape  parameters  of  the 

limbs,  and  h'  will  be  discussed  later.  The  linear  shift  of  C  with  the 

coefficient  of  skew,  SK,  was  calibrated  by  setting  C  =  2  for  zero  skew 

and  C  =  3.5  for  a  coefficient  of  skew  of  4.  Although  this  calibration 

is  judgmental,  it  is  based  on  the  geometrical  restrictions  of  the 

v-scale  in  Fig.  1  and  has,  in  our  experience,  performed  well.  With 

this  calibration,  C  is  defined  by  equation  [2] . 

C  =  2  +  0.375  SK  [2] 


The  boundary  of  v,  VB,  shifts  the  same  amount  to  the  right  as  C  and 
starts  from  VB  =  4.05  for  zero  skew.  We  have  found  this  to  be  a 
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rational  value  beyond  the  required  finite  value  of  4.0  for  the 
calibration  point.  The  expression  for  VB  is  given  in  equation  [3]. 

VB  =  4.05  +  (C  -  2.0)  [3] 

The  h'  in  equation  [1]  is  a  re-scaled  value  of  the  variate  h. 

The  re-scaling  reported  earlier  by  Thcmas  and  Snyder,  1984c,  utilized 
only  the  mean  and  standard  deviation  of  h.  We  later  found  this  to  be 
insufficient  re-scaling  of  the  variate  for  highly  skewed  samples 
because  of  insufficient  reduction  in  curvil inear ity  of  the  sample  when 
plotted  against  probability.  The  solution,  we  found,  was  to  shift  the 
h  and  h'  scales  relative  to  each  other  as  sample  skew  changed.  This 
shifting  is  represented  in  Fig.  1  and  equation  [4] . 

_  h-h  +  a  +  bSK  [4] 

®h 

where  h  is  the  mean  of  the  sample  under  analysis  and  is  its 
standard  deviation.  Equation  [4]  includes  important  elements  of  the 
earlier  version  but  adds  a  shift  controlled  by  skew.  The  a  and  b  are 
calibration  parameters  and  have  been  judgmentally  evaluated  by  making 
the  shift  zero  for  zero  skew  and  one-half  the  standard  deviation  for 
SK  =  4.  Equation  [5]  results  and  may  be  considered  the  operating 
equation  for  re-scaling  h  to  h ' . 

h,  _  h  -  h  |  SK  [5] 

®h  8 

The  last  step  is  to  calibrate  the  shape  parameters,  f  and  d,  of 
the  two  exponential  limbs  given  in  equation  [1] .  This  calibration 
process  requires  the  selection  of  an  hm^  value  which  is  one 
coordinate  of  a  calibrating  point  for  the  v^  limb  and  is  discussed 
below.  The  second  coordinate  for  this  point  is  v^  =  4. 
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Substituting  these  coordinate  values  in  the  limb  of  equation  [1] 

yields  parameter  f  as  shewn  in  equation  [6] . 

f  =  In  ((4  -  VB)/(C  -  VB))/h'  .  [6] 

min 

h'min  is  ccrnputed  by  setting  h  =  h  .  in  equation  [5]  . 

Differentiating  equation  [1]  and  setting  the  derivatives  equal  at  C 
where  h'  =0,  yields  the  solution  for  d  given  in  equation  [7] . 

d  =  f  (VB  -  C)/C  [7] 

The  compound  curve  for  positively  skewed  samples  is  thus  calibrated 
and  its  shape  is  based  on  sample  mean,  standard  deviation,  and 
coefficient  of  skew.  The  only  decision  by  the  user  is  the  selection 
of  a  value  for  h  .  .  To  select  an  appropriate  value,  one  should 
clearly  visualize  the  operating  effect  and  function  of  h^^.  The 
h^  value  sets  the  lower  boundary  of  the  class  containing  the 
smallest  events.  This  corresponds  to  the  requirement  of  setting  the 
first  class  boundary  when  preparing  sample  histograms  in  simplistic 
data  reduction.  The  specification  of  h  .  does  not  establish  a 
lower  limit  of  the  variate  h. 

Negative  Skew 

The  transform  for  negatively  skewed  samples  is  identical  but 

reversed  to  that  represented  in  Fig.  1  for  positive  skew.  Small 

values  of  the  historical  variate,  h,  can  now  extend  to  negative 

infinity.  The  ccmon  point,  C,  moves  to  the  left  with  increasing 

negative  skew.  The  limiting  boundary,  VB,  shifts  from  upper  right  to 

lower  left  in  the  lower  portion  of  Fig.  1.  The  calibration  point 

(v,  =  4,  h  .  )  for  positive  skew  changes  to  (v_  =  0,  h  )  for 
1  nun  ^  ^  2  '  max 

negative  skew  where  h  is  now  an  input.  With  these  differences, 

^  max 

the  equivalent  of  equation  [1]  becomes  equation  [8]  for  negative  skew. 


I 

I 

I 

I 

I 

i 


I 


[ 

I 


7 


V1  =  4 


(4  -  C)  exp  (f  h') ,  v  —  C 
v2  =  (C  +  VB)  exp  (-d  h')  -  VB,  v  <  C 


[8] 


Equation  [2]  which  shifts  the  ccrmon  point,  C,  does  not  change  except 
SK  is  negative.  The  new  equation  for  VB  becomes  equation  [9]. 

VB  =  0.05  +  (2.0  -  C)  [9] 


The  re-scaling  of  h  to  h'  is  unchanged  from  that  shown  in  equation  [5] 
except  for  the  negative  value  of  SK.  Parameter  d  of  equation  [8]  is 
obtained  through  calibration  to  give  equation  [10] . 

d  =  -  In  (VB/(C  +  VB))/h'  [10] 

IllQA 

h'  in  equation  [10]  is  calculated  from  the  input  h  using 
equation  [5] .  Requiring  mathematical  continuity  at  point  C  gives 
parameter  f  as  shewn  in  equation  [11] . 

f  =  d  (C  +  VB)/(4  -  C)  [11] 


A  similar  analogy  exists  for  the  selection  of  h  value  to  that 

max 

given  for  the  selection  of  h  .  . 
r  min 

SMOOTHING 

Once  the  data  have  been  transformed  by  the  above  method,  any 
number  of  smoothing  procedures  could  be  readily  adapted.  However,  we 
oriented  the  transform  development  towards  a  ccmputer-based  form-free 
optimization-smoothing  procedure  known  as  sliding  polynomials.  We 
smooth  the  data  in  the  cumulative  probability  domain,  rather  than  the 
probability  density  domain,  as  demonstrated  in  Fig.  1.  The  smoothing 
procedure  in  this  report  is  similar  to  that  described  by  Thcmas  and 
Snyder,  1984c,  with  a  few  noteworthy  exceptions. 

The  transform  allows  for  reflection  of  any  values  of  the  variate, 
h,  to  values  of  the  intermediate  variate,  v.  With  this 
transformation,  it  is  possible  to  tally  a  sample  of  items,  h,  into 
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classes  which  are  linear  in  v-scale  as  follows.  The  sample  space  in 
v,  0  <  v  -  4  is  divided  into  40  equal  classes,  the  class  boundaries 
are  reflected  to  h-scale ,  and  the  sample  h  is  tallied  into  these  40 
classes.  The  class  subtotals  are  accumulated  frcm  largest  to  smallest 
h.  The  class  accumulations  are  divided  by  the  sample  size,  producing 
accumulated  class  totals  as  ratios  running  from  zero  to  one  across 
four  polynomial  spans,  0-1,  1-2,  2-3  and  3-4.  The  earlier  method 
(Thomas  and  Snyder,  1984c)  used  only  30  classes  with  three  polynomial 
spans.  The  increase  in  the  number  of  classes  and  spans  improved  the 
ability  to  smooth  complex  curvi linearity  found  in  samples  reported  by 
Snyder  and  Thomas,  1986;  and  Thomas  and  Snyder,  1986a,  1986b. 

Since  the  accumulation  is  from  large  h  to  small  h,  the  class 
subtotals  of  probability  must  be  plotted  at  the  right-hand  class 
boundaries  in  v-scale.  This  procedure  obviates  the  long-standing 
problem  of  "plotting  positions"  for  individual  items  of  the  sample. 

Not  all  classes  of  the  40  specified  will  necessarily  be  occupied, 
particularly  for  small  samples.  Consequently,  "plateaus"  occur  in  the 
sample  accumulated  probabilities  which  must  be  smoothed  with  the  arcs 
of  the  sliding  polynomials.  We  should  note  that  smoothing  results 
from  minimizing  errors  in  the  probability  dimension  which  is 
represented  by  the  ordinate  in  the  top  part  of  Fig.  1.  Optimization 
does  not  occur  in  the  scale  of  the  transformed  variate,  v,  which  is 
represented  by  the  abscissa. 
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DISCUSSION 

Two  historical  samples  were  selected  to  demonstrate  the  results 

of  the  transformation  and  smoothing.  The  first  sample  was  a  99-year 

record  of  total  rainfall  during  November  at  Athens,  GA.  The  second 

sample  was  a  47-year  record  of  minimum  temperature  for  November  at  the 

Southern  Piedmont  Conservation  Research  Center  at  Watkinsville ,  GA. 

Table  1  presents  the  sample  characteristics  of  mean,  standard 

deviation,  skew,  and  range  for  the  two  samples.  The  transformation 

parameters  of  C,  VB,  f  and  d  are  also  included.  The  selected  input 

values  of  h  .  and  h  for  the  positive  and  negative  skew, 
min  max  ^  ^ 

respectively,  are  also  shewn. 

The  rainfall  and  temperature  samples  are  plotted  as  accumulated 
probabilities  in  Figs.  2  and  3,  respectively.  Specifically  for  each 
figure,  40  class  ratios  are  plotted  against  the  class  limit  in 
v-scale.  For  the  positive-skew  rainfall  sample,  Fig.  2,  the  smoothing 
line  must  pass  through  v  =  0,  P  =  0  with  zero  slope.  In  contrast  for 
the  negative-skew  temperature  sample.  Fig.  3,  the  smoothing  line  must 
pass  through  v  =  4,  P  =  1  with  zero  slope.  Super- imposed  on  the 
sample  plots  and  derived  smoothing  lines  are  the  v-h  transforms  for 
each  sample.  Those  plotted  transforms  are  derived  from  equations  [1] 
and  [5]  for  positive  skew  in  Fig.  2  and  from  equations  [8]  and  [5]  for 
negative  skew  in  Fig.  3.  The  utility  of  the  v-h  transform  is 
exemplified  by  transformation  of  the  median  value.  Fran  0.5 
probability,  one  goes  horizontally  to  the  smoothing  curve  on  the 
sample,  then  vertically  to  the  v-h  transform,  then  again  horizontally 
to  the  scale  of  the  original  variate  h.  The  mean  of  h  is  also  shewn 
for  canparison. 
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An  APPENDIX  is  included  which  gives  information  and  instructions 
designed  to  assist  a  user  of  this  methodology.  The  section  on  INPUT 
VARIABLES  provides  the  code  for  the  input  variable  name  with  a  brief 
ccrrment.  Likewise,  the  section  on  OUTPUT  VARIABLES  gives  the  code  for 
the  output  variable  name  with  a  brief  cement.  Next  is  the  program 
listing,  sectionalized  with  remark  statements.  Last,  is  the  sample 
output  for  the  analysis  of  the  November  rainfall  which  had  a  positive 
skew.  Not  shown  is  the  sample  output  for  the  November  temperature. 

In  conclusion,  this  research  report  presents  a  variate 
transformation  that  appears  to  have  great  potential  for  use  with  many 
types  of  stochastic  distributions.  The  method  is  oriented  toward  the 
piece-wise  method  of  sliding  polynomials;  however,  it  is  not  limited 
to  this  particular  methodology.  The  method  has  proved  useful  for 
simple  samples  as  well  as  more  complex  samples  exhibiting  such 
characteristics  as  bi -modal  tendency,  positive  or  negative  skew,  and 
presence  or  absence  of  outliers. 
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TABLE  1.  CHARACTERISTICS  OF  TEST  DATA 


November  Total  November  Kinimum 

Monthly  Rainfall  Temperature 


* 


Mean 

3.95 

(7.77) 

23.87 

(4.52) 

Std.  Dev.* 

2.33 

(5.92) 

5.08 

(2.82) 

Coef.  of  Skew 

2.06 

-0.99 

★ 

Range 

15.78- 

(40.08- 

35- 

(1.67- 

0.33 

0.84) 

5 

-15.00) 

c** 

2.77 

1.63 

VB** 

4.82 

0.42 

f** 

0.86 

0.61 

d** 

0.64 

0.70 

h  .  * 
nun 

0 

(0) 

— 

* 

h 

max 

— 

36.00 

(2.22) 

* 

Units :  rainfall 

,  inch 

(cm)  ;  temperature,  °F 

(°C). 

♦♦Definitions:  C  -  camon  point,  VB  -  asymptotic  boundary,  f  and 
d  -  shape  parameters. 


Fig.  1 


Schematic  for  variate  transformation  and  smoothing. 
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v  -  Scale 


Fig.  2  Transformed  and  smoothed  sample  of  total 

monthly  rainfall  for  November  at  Athens,  GA. 


Rainfa  1 1 ,  c  m 


Probability  (v  <  v 
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Fig.  3  Transformed  and  smoothed  sample  of  minimum 

monthly  temperature  for  Nov  at  Watkinsville,  GA. 
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APPENDIX 


The  Appendix  includes  a  description  of  input  and  output  variables, 
program  listing,  and  sample  output  of  the  analysis  of  the  November  total 
rainfall. ^ 


2 

The  program  is  presented  for  the  convenience  of  potential 
users.  While  the  program  has  been  run  and  tested  on  various  data 
sets,  the  originators  of  the  program  assume  no  responsibility  for 
its  accuracy  or  adequacy.  Such  responsibility  must  rest  solely  on 
the  user.  We  stand  ready  to  assist  and  advise  within  the 
limitations  imposed  by  our  operating  resources.  The  program  is 
listed  in  Hewlett-Packard  BASIC  with  Ext.  2.1.  (Trade  name  is 
included  for  the  benefit  of  the  reader  and  does  not  imply  an 
endorsement  or  preferential  treatment  of  the  named  product.) 
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Input  Variables 

Variable  name 

Cement 

PRTR 

Printer  output  device  control. 
(1=€KT,  706=Printer) 

T$ 

Problem  title 

HM 

Minimum  class  limit  in  data  scale 

HH 

Maximum  class  limit  in  data  scale 

DIC 

Data  input  control 

(l=Keyboard  input,  2=Disk  resident) 

N 

Number  of  data  points 

DSN$ 

Disk  resident  data  file  name 

H (I) ;  1=1, N 

Data  points 

19 


Output  Variables 

Variable  Name 

Canments 

T$ 

Problem  title 

HM 

Minimum  class  limit  in  data  scale 

HH 

Maximum  class  limit  in  data  scale 

HB 

Sample  mean 

SSD 

Sample  standard  deviation 

SK 

Sample  coefficient  of  skew 

Cl 

Carmen  point  of  exponential  limbs 

VB 

Asymptotic  boundary 

F 

Shape  parameter 

D 

Shape  parameter 

CH (I) : 1=1 ,40 

Class  limit  in  h-scale 

CT(I) ; 1=1, 40 

Sample  class  probabilities 

SC (1,6) : 1=1 , 5 

'Best  Fit'  nodal  ordinates 

CP (I) : 1=1, 40 

Smoothed  class  probabilities 

SD (I) :I=1,5 

Standard  deviation  of  nodes 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 


20 


Program  Listing 

RIM — - ID  PROBABILITY  ANALYSIS- - 

REM - - 4  SPANS  WITH  NEW  TRANSFORM - - 

REM - —HANDLES  POS  AND  NEG  SKEW - - 


OPTION  RASE  1 

DIM  H(100)  ,CH (40)  ,CT(40)  ,0(40,7)  ,SC(6,6)  ,CP(40)  ,E(40)  ,SD(5) 
DIM  T$ [50] ,DSN$  [30] 

INPUT  "TO  HARDCOPY  PRINT  ENTER  ( 706 )  ,  CRT  ENTER  ( 1)  "  ,PRTR 
LINPUT  "  PROBLEM  TITLE ",T$ 

PRINT  "MIN  H  USED  WITH  0  OR  +  SKEW,  MAX  H  USED  WITH  -  SKEW" 
INPUT  "MINIMUM  VALUE  CLASS  LIMIT  OF  H"  ,HM 
INPUT  "MAXIMUN  VALUE  CLASS  LIMIT  OF  H"  ,HH 

RE>1____INPUT  DATA  CONTROL  IS  '1'  FOR  KEY-IN  '2'  FOR  DISK- 
INPUT  "DATA  CONTROL:  ENTER(l)  FOR  KEY-IN,  (2)  FOR  DISK",DIC 
OUTPUT  PRTR; "  NEW  TRANSFORM  WITH  FOUR  ARCS" 

OUTPUT  PRTR;T$ 

OUTPUT  PRTR;  "MINIMUM  VALUE  CLASS  LIMIT  OF  H  ";HM 

OUTPUT  PRTR;  "MAXIMUN  VALUE  CLASS  LIMIT  OF  H  "  ;HH 

IF  DIC=2  THEN  240 

INPUT  "NUMBER  OF  DATA  POINTS ",N 

FOR  1=1  TO  N 

INPUT  H (I) 

NEXT  I 
GOTO  310 

LINPUT  "  DATA  SET  NAME  ",DSN$ 

ASSIGN  @INP  TO  DSN$ 


ENTER  @INP;N 


270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 
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FOR  1=1  TO  N 
ENTER  @INP;H (I) 

NEXT  I 

ASSIGN  @INP  TO  * 

REM - — • - FIND  FIRST  THREE  MOMENTS - - 

SH=0. 

SHH=0 . 

SHHH=0 . 

FOR  1=1  TO  N 
SH=SH+H ( I ) 

XH=H (I) *H (I) 

SHH=SHH+XH 
SHHH=SHHH+XH*H (I) 

NEXT  I 
HB=SH/N 

SSD=SQR (SHH/N-HB*HB) 

SK= ( (SHHH/N) - (3*HB*SHH/N) +(2*HB*HB*HB) ) /SSD/SSD/SSD 

OUTPUT  PRTR; "  MEAN  "  ;HB; "  SID  DEV  ";SSD;"  COEF  OF  SKEW  ";SK 


REM - TEST  FOR  NEGATIVE  OR  POSITIVE  SKEW - - - 

IF  SK<0.  THEN  810 

REM - POSITIVE  SKEW - 


HPM=(HM-HB) /SSD+SK/8 
C1=2.+.375*SK 
VB=4 .05+ (Cl-2 . ) 

F=LOG ( (4-VB) / (Cl-VB) ) /HFM 

D=F* (VB-C1) /Cl 

V=0. 


540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 
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OUTPUT  PRTR;"C";C1;"  VB";VB;"  F";F;"  D";D 

REM— — - — CALCULATE  CLASS  LIMITS— - - 

FOR  1=1  TO  39 
V=V+.l 

IF  V>C1  THEN  610 

CH (I) =-SSD* (LOG (V/Cl) /D+SK/8) +HB 
GOTO  620 

CH (I) =SSD* (LOG ( (V-VB) / (Cl-VB) ) /F-SK/8) +HB 
NEXT  I 
CH  (40)  =HM 

REM- — - - - - -CALCULATE  SAMPLE  PROBABILITIES' 

FOR  J=1  TO  N 

FOR  1=1  TO  40 

IF  H (J) >=CH (I)  THEN  690 

GOTO  710 

CT  (I)  =CT  (I)  +1 

GOTO  720 

NEXT  I 

NEXT  J 

FOR  J=2  TO  40 
CT(J)=CT(J)+CT(J-1) 

NEXT  J 

FOR  1=1  TO  40 
CT  (I)  =CT  (I)  /N 
NEXT  I 
GOTO  1140 


REM- 


•NEGATIVE  SKEW- 


810  HPM=(HH-HB)/SSD+SK/8 
820  C1=2.+.375*SK 

830  VB=.05+(2.-Cl) 

840  D=-LDG (VB/ (Cl+VB) ) /HPM 
850  F=D* (Cl+VB) / (4-C1) 

860  OUTPUT  PRTR; "C" ;C1 ; "  VB";VB;"  F";F;"  D";D 

870  RIM- - - - - - CALCULATE  CLASS  LIMITS- - - 

880  V=0 . 

890  FOR  1=2  TO  40 

900  V=V+.l 

910  IF  V>C1  THEN  940 

920  CH  (I)  =SSD*  (-LOG  ( (V+VB)  /  (Cl+VB) )  /D-SK/8)  +HB 
930  GOTO  950 

940  CH  (I)  =SSD*  (LOG  ( (4-V)  /  (4-C1) )  /F-SK/8)  +HB 
950  NEXT  I 
960  CH (1) =HH 

970  RIM - CALCULATE  SAMPLE  PROBABILITIES  — — 

980  FOR  J=1  TO  N 

990  FOR  1=2  TO  40 

1000  IF  H (J) >=CH (I)  THEN  1020 

1010  GOTO  1040 

1020  CT(I)=CT(I)+1 

1030  GOTO  1050 

1040  NEXT  I 

1050  NEXT  J 

1060  FOR  J=2  TO  40 

1070  CT (J)  =CT ( J) +CT (J-l) 


1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 
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NEXT  J 

FOR  1=2  TO  40 
CT(I)=CT(I)/N 
NEXT  I 
CT  (1) =0 . 

OUTPUT  PRTR;  "CLASS  LIMITS  AND  SAMPLE  CLASS  PROBABILITIES" 
J=1 

FOR  1=1  TO  40 

IMAGE  #,3D,7D.DD,2D.3D 

OUTPUT  PRTR  USING  1170;I,CH(I) ,CT(I) 

IF  J=3  THEN  1230 
J=J+1 

OUTPUT  PRTR;" 

GOTO  1250 
J=1 

OUTPUT  PRTR 
NEXT  I 
OUTPUT  PRTR 

REM - LOCATE  SEVEN  NODES— - - - - — — 

XN  (1)  =1 
FOR  1=2  TO  7 
XN (I) =XN (1-1) +1 
NEXT  I 


RZM - CALCULATE  COEFFICIENT  C(I,J) 

REM - 7  FOR  EACH  OF  40  CLASSES- 


FOR  IP=1  TO  40 


1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 
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JC=1 

IF  (IP*.1)<=JC  THEN  1390 
JC=JC+1 
GOTO  1360 
DP=(IP*.l) -JC+1. 

Z=- . 5+DP 

IF  SK<0.  THEN  Z=Z-.l 
C (IP, JC) =( ( (-8*Z+4) *Z+2) *Z-1) /16 
C(IP,JC+1)  =  (  ( (24*Z-4)*Z-22) *Z+9) /16 
C  (IP, JC+2)  =  ( ( (-24*Z-4) *Z+22) *Z+9)  /16 
C  (IP, JC+3)  =  ( ( (8*Z+4) *Z-2)  *Z-1)  / 16 
NEXT  IP 

IF  SK<0.  THEN  1620 

HEM - IF  POS.  SKEW,  PUT  IN  LEFT  BOUNDARY  ZERO  AND— 

REM - -ZERO  SLOPE,  RIGHT  BOUNDARY  IS  FREE - — - - — 

FOR  IP=1  TO  40 
C(IP,3)=C(IP,3)+C(IP,1) 

NEXT  IP 

REM- - MOVE  TO  LEFT  AND  ADD  CLASS  PROB  TO  MATRIX- — - - 

FOR  IP=1  TO  40 
FOR  J=1  TO  5 
C (IP,  J)  =C  (IP, J+2) 

NEXT  J 

C  (IP,6)  =CT(IP) 

NEXT  IP 
GOTO  1720 

REM— — !F  NEG.  SKEW,  POT  IN  RIGHT  BOUNDARY  1.0  AND - 


■ZERO  SLOPE,  LEFT  BOUNDARY  IS  FREE - 


1620  REM- 
1630  FOR  IP=1  TO  40 


1640  C(IP,5)=C(IP,5)-K:(IP,7) 

1650  NEXT  IP 

1660  REM— ADD  CLASS  PROBABILITIES  AND  NODE  6=1.0 - 

1670  REM - BOUNDARY;  SAVE  C(IP,6)  IN  7 - 


1680  FOR  IP=1  TO  40 
1690  C (IP, 7) =C (IP, 6) 

1700  C  (IP, 6)  =CT  (IP) -C (IP, 7) 

1710  NEXT  IP 

1720  REM - - CALCULATE  SUMS  OF  SQUARES  MATRIX- 

1730  FOR  K=1  TO  40 

1740  FOR  1=1  TO  6 

1750  FOR  J=I  TO  6 

1760  SC(I,J)=SC(I,J)4C(K,I)*C(K,J) 

1770  NEXT  J 
1780  NEXT  I 
1790  NEXT  K 

1800  PRINT  "SUMS  OF  SQUARES  CALCULATED" 

1810  REM - - FILL  IN  COMPLETE  MATRIX - 

1820  FOR  1=1  TO  5 
1830  FOR  J=1  TO  5 
1840  SC(J,I)=SC(I,J) 

1850  NEXT  J 
1860  NEXT  I 

1870  REM - 

1880  FOR  J=1  TO  5 


— SOLVE  SET  OT  5  EQUATIONS' 
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1890  FOR  W  TO  5 

1900  IF  SC(I,J)<>0.  THEN  1940 

1910  NEXT  I 

1920  PRINT  "NO  UNIQUE  SOLUTION" 

1930  GOTO  2180 
1940  FOR  K=1  TO  6 
1950  X=SC (J,K) 

1960  SC(J,K)=SC(I,K) 

1970  SC  (I  ,K)  =X 
1980  NEXT  K 
1990  Y=1/SC (J, J) 

2000  FOR  K=1  TO  6 
2010  SC ( J,K) =Y*SC (J/K) 

2020  NEXT  K 
2030  FOR  1=1  TO  5 
2040  IF  I=J  THEN  2090 
2050  Y=-SC(I,J) 

2060  FOR  K=1  TO  6 

2070  SC(I,K)=SC(I,K)+Y*SC(J,K) 

2080  NEXT  K 
2090  NEXT  I 
2100  NEXT  J 

2110  OUTPUT  PRTR; "  SOLUTION  FOR  FIVE  BEST  FIT  NODES" 

2120  FOR  1=1  TO  5 

2130  OUTPUT  PRTR  USING  "3D,K,DD. 5D,K" ;I , "  ",SC(I,6),"  " 

2140  NEXT  I 
2150  OUTPUT  PRTR 


2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 
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OUTPUT  PRTR 

REM - LAY  IN  SLIDING  POLYNOMIAL  THROUGH  NODES - 

OUTPUT  PRTR; "SMOOTHED  CLASS  TOTAL  PROBABILITIES " 

IF  SK<0.  THEN  2380 
FOR  1=1  TO  40 
IF  I>30  THEN  2330 
IF  I>20  THEN  20 
IF  I>10  THEN  2260 

CP(I)=C(I,1)*SC(1,6)4C(I,2)*SC(2,6) 

GOTO  2360 

CP  (I)  =C  (1 , 1)  *SC  (1  ,  6)  4C  (1 , 2)  *SC  (2, 6)  4C  (1 , 3)  *SC  (3 ,6) 

GOTO  2360 
CP(I)=0. 

FOR  J=1  TO  4 

CP(I)=CP(I)+C(I,J)  *SC (J,6) 

NEXT  J 
GOTO  2360 
FOR  J=2  TO  5 

CP(I)=CP(I)4C(I,J)  *SC  (J,6) 

NEXT  J 

NEXT  I 

GOTO  2500 

FOR  1=1  TO  40 

IF  I>30  THEN  2480 

IF  I>20  THEN  2460 

IF  I>10  THEN  2440 

CP  (I)  =€  (1 , 1)  *SC  (1 ,6)  +C  (1,2)  *SC  (2,6)  +C  (1 ,3)  *SC  (3 ,6) +C (I ,4)  *SC  (4,6) 


2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 
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GOTO  2490 

CP  (I)  (1, 2)  *SC  (2,6) +C (1,3)  *SC  (3 ,6)  +C  (1,4)  *SC  (4,6)  +C  (1,5)  *SC  (5 ,6) 

GOTO  2490 

CP  (I)  ^  (1,3)  *SC  (3 ,6)  +C  (1 ,4)  *SC  (4 ,6)  +C  (1,5)  *SC  (5,6)  -K3  (1,7) 

GOTO  2490 

CP (I)  =C (1,4)  *SC  (4 ,6)  +C (1,5)  *SC (5,6)  +C (1, 7) 

NEXT  I 
J=1 

FOR  1=1  TO  40 

OUTPUT  PRTR  USING  "#,3D,3D.3D,2X";I,CP(I) 

IF  J=5  THEN  2570 
OUTPUT  PRTR;" 

J=J+1 
GOTO  2590 
J=1 

OUTPUT  PRTR 
NEXT  I 

FOR  1=1  TO  40 
E  (I)  =CP  (I) -CT  (I) 

NEXT  I 

FOR  J=1  TO  5 

Z2S=0. 

Z1S=0. 

Z0S=0. 

FOR  1=1  TO  40 

Z2S=Z2S+E (I) *E (I) *C (I , J) *C(I,J) 

zis=zis-k:(i,j)*c(i,j) 
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2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 


Z0S=Z0S+C(I,J) 

NEXT  I 

IF  Z0S<=0.  THEM  2740 
GOTO  2760 
SD (J) =99.99 
GOTO  2770 

SD ( J) =SQR (Z2S/Z1S/Z0S*N/ (N-5) ) 

NEXT  J 

OUTPUT  PRTR; "  STANDARD  DEVIATION  OF  NODES" 
FOR  J=1  TO  5 

OUTPUT  PRTR  USING  "3D,3D.4D";J,SD(J) 

NEXT  J 
OUTPUT  PRTR 

OUTPUT  PRTR;"  PROBLEM  END  " 

STOP 

END 
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Program  Output 


NEW  TRANSFORM  WITH  FOUR  ARCS 
NOV  TOTAL  RAIN 

MINIMUM  VALUE  CLASS  LIMIT  OF  H  0 
MAXIMUN  VALUE  CLASS  LIMIT  OF  H  0 

MEAN  3.06353535354  STD  DEV  2.32778010743  COEF  OF  SKEW  2.0566512138 
C  2.77124420517  VB  4.82124420517  F  .863814250258  D  .63899789478 
CLASS  LIMITS  AND  SAMPLE  CLASS  PROBABILITIES 


1 

14.57 

.010 

2 

12.04 

.010 

3 

10.56 

.010 

4 

9.52 

.010 

5 

8.70 

.010 

6 

8.04 

.040 

7 

7.48 

.040 

8 

6.99 

.051 

9 

6.56 

.071 

10 

6.18 

.091 

11 

5.83 

.121 

12 

5.51 

.141 

13 

5.22 

.152 

14 

4.95 

.162 

15 

4.70 

.172 

16 

4.47 

.182 

17 

4.25 

.232 

18 

4.04 

.253 

19 

3.84 

.273 

20 

3.65 

.333 

21 

3.48 

.374 

22 

3.31 

.394 

23 

3.14 

.404 

24 

2.99 

.404 

25 

2.84 

.404 

26 

2.70 

.444 

27 

2.56 

.495 

28 

2.43 

.495 

29 

2.29 

.525 

30 

2.15 

.556 

31 

1.99 

.606 

32 

1.83 

.646 

33 

1.66 

.727 

34 

1.48 

.778 

35 

1.28 

.828 

36 

1.07 

.869 

37 

40 

.84 

0.00 

.879 

1.000 

38 

.59 

.919 

39 

.31 

.970 

SOLUTION  FOR  FIVE  BEST  FIT  NODES 

1  .08569 

2  .31666 

3  .56916 

4  .98832 

5  .75764 


SMOOTHED  CLASS  PROBABILITIES 


1 

.001 

2 

.004 

3 

.009 

4 

.015 

5 

.023 

6 

.033 

7 

.044 

8 

.057 

9 

.070 

10 

.086 

11 

.103 

12 

.122 

13 

.144 

14 

.167 

15 

.191 

16 

.216 

17 

.241 

18 

.267 

19 

.292 

20 

.317 

21 

.340 

22 

.363 

23 

.386 

24 

.408 

25 

.431 

26 

.455 

27 

.480 

28 

.508 

29 

.537 

30 

.569 

31 

.607 

32 

.653 

33 

.703 

34 

.756 

35 

.809 

36 

.859 

37 

.905 

38 

.943 

39 

.972 

40 

.988 

STANDARD  DEVIATION  OF  NODES 

1  .0040 

2  .0071 

3  .0048 

4  .0076 

5  99.9900* 

PROBLEM  END 


♦Default  for  undefined  value 
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